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This  paper  presents  a  simple,  analytical  expression  for  A  ,  derived  by  assuming  the 
lower  truncation  limit  to  be  zero.  With  an  independent  estimate  of  the  value  of  the  lower 
truncation  limit ,  a  factor  f  is  also  determined ,  which  will  correct  A  to  99  percent’  pf  its 
actual  value. 
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The  Application  of  Doubly-Truncated 
Hydrometeor  Distributions  to  Numerical 

Cloud  Models 

1.  INTRODUCTION 

The  presence  of  cloud  or  precipitation  particles  or  "hydrometeors"  (both  liquid  and 
frozen)  has  a  significant  impact  on  many  Air  Force  systems  and  operations.  Because  of 
this,  it  is  often  necessary  to  determine  reasonable  particle  size  distributions  from  limited 
data,  which  may  include  measurements  of  mean  radar  reflectivity  or  mean  densities  of 
liquid  water  or  ice  over  a  region  of  the  atmosphere.  An  important  tool  for  studying  the 
evolution  of  cloud  properties  and  processes  is  the  numerical  cloud  model.  When  clouds 
are  modeled  in  two  or  three  dimensions,  it  is  again  often  necessary  to  obtain  size  distri¬ 
butions  of  liquid  water  or  ice  particles  given  only  data  averaged  over  a  region  of  space 
or  "grid  volume".  The  averaged  dsta  generally  include  the  mixing  ratios  of  the  various 
hydrometeor  species  (cloud  water,  rain,  cloud  ice,  snow  or  aggregates,  graupel/hail, 
etc.)  from  which  one  can  calculate  the  density  ("liquid  water  content"  or  "ice  wster  con¬ 
tent")  of  each  species.  In  the  cloud  model  it  is  necessary  to  know  the  distribution  of 
each  particle  type  in  order  to  derive  mathematical  expressions  for  the  conversion  pro¬ 
cesses  among  the  various  species  (for  example  the  rate  of  freezing  of  raindrops  to  form 
graupel  depends  on  the  size  distribution  of  the  raindrops).  Once  a  suitable  relationship 
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between  mean  cloud  properties  and  size  distribution  is  derived  and  validated,  it  is  also 
useful  for  including  into  diagnostic  cloud  modeling  and  artificial  intelligence  schemes 
involving  clouds. 

The  Cloud  Physics  Branch  of  the  Atmospheric  Sciences  Division  at  the  Air  Force 

Geophysics  Laboratory  (AFGL)  has  been  utilizing  a  3- dimensional,  numerical-cloud 

model  (hereafter  referred  to  as  the  3-D  model)  for  predicting  mixing  ratios  of  hydro- 

meteor  liquid  and  ice  water  content  from  bulk  microphysical  parameterizations.  This 

1  2 

AFGL  3-D  model  is  a  version  of  that  developed  at  Colorado  State  University.  ’  The 

parameterizations  rely  upon  assumed  exponential  size  distributions  for  deriving  the 

conversion  processes  occurring  amongst  the  various  hydrometeor  species. 

The  assumption  of  exponential  distributions ,  where  the  logarithmic  values  of  the 

number  of  hydrometeors  decrease  as  a  function  of  increasing  diameter  size,  is  generally 

3-7 

accepted  as  being  reasonably  descriptive  and  is  supported  by  various  studies.  Fig¬ 
ure  la  shows  a  typical  semi- logarithmic,  number-density  plot  where  Figure  lb,  the  same 
data  plotted  linearly,  emphasizes  the  magnitude  differences  in  numbers  of  hydrometeors 
at  the  small  and  large  diameters. 

In  the  current  modeling,  the  assumed  exponential  distributions  are  generally  de¬ 
fined  by  constant  values  of  the  total  number  of  particles  or  for  the  slope  and  /or  zero 
intercept  of  the  semi-logarithmic  number  density  plot,  although  observations  do  not  sup¬ 
port  these  generalities.  Derived  distributions  employing  these  types  of  assumptions 
have  very  little  chance  of  agreeing  with  measurements  of  actual  events. 
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_  .aa, 

3.  Marshall,  J.S.,  and  Palmer,  W.  McK.  (1948)  The  distribution  of  raindrops  with 

size,  J.  Meteorol.  5:165-166. 

v 

4.  Marshall,  J.S.,  and  Gunn,  K.L.S.  (1952)  Measurement  of  snow  parameters  by 

radar,  J,  Meteorol.  8:322. 

vV 

5.  Iraai,  I.,  Fujiwara,  M.,  Ichimura,  I.,  and  Toyama,  Y.  (1955)  Radar  reflectivity 

of  falling  snow,  Pap,  in  Meteorol.  and  Geophys.  (Japan)  6:130-139. 

6.  Gunn,  K.L.S. ,  and  Marshall,  J.S.  (1958)  The  distribution  with  size  of  aggregate 

snowflakes ,  J.  Meteorol.  1^: 452  ( 47 9) . 


7.  Ohtake,  T.,  and  Henmi,  T.  (1970)  Radar  reflectivity  of  aggregated  snowflakes, 
Preprints  of  papers  presented  at  the  14th  Radar  Meteorology  Conference. 
Tucson,  Arizona,  17-20  November  19^0,  pp  209-211. 


Diameter 


Diameter 


Figure  la.  A  Typical  Exponential 
Number  Distribution  Plotted  in  a 
Semi-Logarithmic  Format 


Figure  lb.  A  Typical  Exponential 
Number  Distribution  Plotted  in  a 
Linear  Format 


Constant  values  for  the  slopes  of  exponential  distribution  lines  are  not  realistic 
assumptions,  as  considerable  variations  in  slope  values  are  evident  even  in  measurements 
of  the  same  hydrometeor  type.  One  particular  rain  situation,  for  example,  has  shown 

slope  values  that  range  between  1.06  and  2.2  mm  1  for  distributions  averaged  over  a 
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r-l  m  volume.  More  pronounced  differences  in  number  density  slopes  occur  with 
changes  of  hydrometeor  type  (rain,  snow,  ice  crystals,  etc.)  with  rain  generally  having 
the  smallest  values  and  the  high-altitude  ice  crystals  the  largest.  A  considerable  over¬ 
lapping  of  values  is  apparent  between  the  various  types  and  mixtures  of  types,  prevent¬ 
ing  specific  distribution  slopes  being  assigned  to  particular  hydrometeors. 

The  study  mentioned  above  also  shows  number  totals  of  precipitable  sized  drops 
ranging  from  227  to  1239  for  the  same  averaging  period.  Measurements  such  as  these 
point  to  substantial  variations  in  both  the  slopes  of  the  distribution  lines  and  the  zero 
intercept  of  the  number  density  plots,  as  both  these  parameters  have  an  effect  on  the 
total  number  of  hydrometeors. 

The  amount  of  water  contained  in  a  distribution  and  the  factor  of  radar  microwave 
energy  reflected  from  the  distributed  hydrometeors  may  be  calculated  from  an  assumed 
size  distribution.  Figure  2a  shows  the  distributed  water  content  and  Figure  2b,  the 


8.  Berthel,  R.O.,  and  Plank,  V.G.  (1983)  A  model  for  the  estimation  of  rain  distri¬ 
butions,  Environmental  Research  Papers,  No.  822,  AFGL-TR-83-0030. 
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Figure  2a.  A  Plot  of  the  Distrib¬ 
uted  Liquid  Water  Content  from 
the  Number  Density  Distribution 
of  Figure  1 


Figure  2b.  A  Plot  of  the  Distrib¬ 
uted  Radar  Reflectivity  Factor 
from  the  Number  Density  Distri¬ 
bution  of  Figure  1 


distributed  radar  reflectivity  factor  that  would  result  from  a  full  exponential  distribu¬ 
tion  of  rain  drops  as  in  Figure  1.  The  plots  indicate  that  not  only  do  these  parameters 
peak  at  different  drop  sizes  but  that  the  relative  water  contribution  (amount  at  a  spe¬ 
cific  diameter  as  compared  to  the  total)  is  greater  than  that  from  the  relative  radar 
reflectivity  factor  at  the  smaller  sizes ,  and  the  opposite  effect  is  prevalent  at  the  larger 
diameters.  Both  parameters  show  minor  relative  contributions  at  the  very  small  diam¬ 
eters  with  more  substantial  increases  from  the  lesser  numbers  of  larger-sized  drops. 

In  this  report,  an  exponential  distribution  refers  primarily  to  precipitation -sized 
hydrometeors,  as  the  small  cloud  and  aerosol  particles  may  form  separate  distributions 
of  varying  shapes  and  numerical  concentrations  depending  upon  atmospheric  conditions 
(Section  5).  The  separation  of  cloud  and  precipitation  hydrometeors  is  normally  con¬ 
sidered  to  occur  between  .05  and  .1  mm  in  drop  diameters.  Because  of  the  exponential 
increase  in  number  with  decreasing  diameter,  the  total  number  (per  unit  volume)  of 
precipitable  particles  is  highly  dependent  upon  this  lower  truncation  limit ,  although  the 
sensitivity  of  total  water  content  or  radar  reflectivity  to  this  parameter  is  small.  On  the 
other  hand,  truncating  or  terminating  exponential  distributions  at  the  larger  diameters 
has  very  little  effect  on  total  number  but  makes  considerable  differences  in  the  total 
water  content  and  radar  reflectivity  factor. 


'.vvV-iJ 
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The  following  sections  describe  s  doubly  truncated,  hydrometeor  model  of  ex¬ 
ponential  type  that  allows  variability  in  all  the  distribution  parameters  with  only  the 
shape  of  the  number  density  curve  being  assumed .  This  hydrometeor  model  is  then 


combined  with  empirically  derived  relationships  to  tailor  distributions  to  particular 
atmospheric  situations. 

2.  EXPONENTIAL  DISTRIBUTION  EQUATIONS 

It  has  been  demonstrated  by  several  authors  that  the  size  distribution  proper¬ 
ties  of  raindrops,  snowflakes,  and  ice  crystals  can  be  reasonably  described  by  a 
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distribution  function  of  exponential  type.  This  distribution  function  specifies 
that  the  number  concentration  of  the  hydrometeor  particles  will  decrease  with  increas¬ 
ing  diameter  (or  equivalent-melted  diameter)  in  the  manner 


N  =  NQ  e~°A  No.  m  3  mm  1 


(1) 


where  Nq  is  the  number  per  cubic  meter  per  millimeter  bandwidth  at  the  zero  intercept 

of  the  semi- logarithmic,  number-density  plot,  A  is  the  slope  of  the  exponential  number- 

density  distribution ,  and  D  is  the  drop  diameter  in  millimeters . 

In  the  present  study  Eq .  (1)  is  presumed  to  apply  only  between  the  truncation 

limits  D  =  d  (a  minimum  diameter)  and  D  =  D_  (a  maximum  diameter) .  The  double 

m 

truncation  of  an  exponential  distribution  function  has  been  previously  discussed  by 

g 

Sekhon  and  Srivastava. 

2.1  Number  Total 

The  total  number  of  hydrometeors  (NT)  in  a  population  described  by  Eq.  (1)  is 


N_  = 


dD  No.  m 


-3 


(2) 


or 


N  =  to 

nt  a 


No.  m 


-3 


(3) 


Sekhon,  R.S.,  and  Srivastava,  R.C.  (1970)  Snow  size  spectra  and  radar  reflec¬ 
tivity.  J.  Atmos.  Sci.  27:299-307. 

-  v\ 


The  liquid-water  content  or  mass  (M)  of  the  hydrometeor  populations  per  unit  vol 
ume  is  distributed  with  diameter  as 


m  -  JLV  1/1 — ^  o  m  n  ^  D  A  -3  -i 

Mjj  —  g  x  10  P  ^  N  Q  D  e  ^  ill  mm 


where  p_  is  the  density  of  liquid  water  in  g  cm 


The  total  liquid-water  content  of  the  population  is 


M  =/m  Mn  dD  g  nf 
d 


that,  from  Eq.  (6)  and  integration,  yields 


*x  1°  p  w  N0  r  (4)rM  -3 

M  =  - - -  g  m 

6  A  4 


where  r  (4)  is  the  gamma  function  of  4  and  r^  is  a  truncation  ratio  for  liquid -water 
content  given  by 


/mMD  dD 
MD  dD 


M 


:-|-"£e_<1A  [(dA)3  +  3(dA)2  +  6  dA  +  6]  - 


e'®mA  [(D„A)3  +  3(DmA)2  +  6DmA  +  6]  J 


m 


m 


(10) 


2.3  Radar  Reflectivity  Factor 

The  distributed  values  of  the  equivalent  radar  reflectivity  factor  (Z)  for  the  hy¬ 
drometeor  populations  are  specified  by 


Zp  =  Ng  D®  e~^A  mm®  m  3  mm  *  . 


(11) 


The  total  value  of  the  radar  reflectivity  factor,  for  the  entire  population  of  hy¬ 
drometeors,  is 


Z  =  /  m  ZD  dD  mm6  m  3 
d 


or,  from  Eq.  (11)  on  integration, 


(12) 


Z  - 


Nq  T  ( 7)  rz  mm6  m  3 


(13) 


where  r  (7)  is  the  gamma  function  of  7  and  r^  is  the  truncation  ratio  for  the  radar  re¬ 
flectivity  factor  as  defined  by 


that  becomes 


=  J_  f 

720  1« 


((dA)6  +  6(dA)5  +  30(dA)4  + 


120(dA)°  +  360(dAr  +  720dA  +  720] 


-e-®mA  I- (DA)6  +  6(DmA)5  +  30(D  A)4  + 
l  m  m  m 


120(D  A)3  +  360(D  A)2  +  720(D  A )  +  720]  |  . 
m  m  m  l 


2.4  Modal  Diameters 


The  "modal  diameters"  of  the  MQ  and  distributions  are  characteristic  param¬ 
eters  of  the  hydrometeor  populations.  These  diameters,  which  specify  the  maximum 
value  points  of  liquid -water-content  and  radar- reflectivity  factor,  are  given  respec¬ 
tively  by 


D'm  =  3/A  mm 


D'^,  =  6/A  mm. 


2.5  Median  Volume  Diameter 


An  additional  characteristic  parameter  of  the  distribution  which  is  conven¬ 


tionally  used  and  referenced  is  the  median  volume  diameter  (Dg).  This  diameter  satis¬ 


fies  the  integral  relation 


Dn  D 

Md  dD  =/m  Md  dD 
d  Dq 


that,  when  the  integration  is  performed  using  Eq.  (6),  and  if  all  A  terms  are  included 
on  the  right,  yields 
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F %r<r ’.r*/  t •  r* *  •  ^ . «r» r, ^ ,•* .  ^ .  *  *  * ,  * ,  • , • , ••  < , 


i  f _  _ 2l(DQA)*  + 

D°  *  L  e-dAKdA)3  +  3(dA)2  +  6dA  +  6]  + 

3(DqA)2  +  6DqA  +  6) _ 

e~®mAf(D  A)3  +  3(DA)2  +  6  D  A  +  6) 
in  m  n 


(19) 


It  is  seen  that  DQ  is  not  separable  in  this  equation.  However,  the  equation  can 
be  readily  solved  by  "trial  and  error"  iteration  once  information  about  d ,  and  A  is 

available. 


2.6  Distribution  Slope 

The  parameter  NQ  may  be  eliminated  between  Eqs.  (8)  and  (13)  to  provide  an  ex¬ 
pression  for  the  "exponential  slope"  of  the  distribution  function  of  Eq .  ( 1)  as 


A  =  61.2 


1/3 


mm 


-1 


(20) 


When  d  and  Dm  are  known,  in  Eqs.  (10)  and  (15),  this  equation  may  also  be 
solved  by  iteration. 

2  7  Coefficient  of  the  Exponential  Distribution  Function 

Nq  may  then  be  determined  through  Eqs.  (8)  or  (13)  as 


N0  = 


M  A 


*x  10  3  r, 


M 


M  -3  -1 

No.  m  mm 


(21) 


or 


N  =-£A_ 
0  720  r„ 


No.  m 


-1 

mm 


(22) 


using  =  1  g  cm-3,  l~(4)  =6  and  F(7)  =  720. 

Equations  (1)  through  (22)  constitute  a  descriptive  equation  set,  or  model,  that 
may  be  used  to  determine  the  parameters  of  hydrometeor  distributions  when  suitable 
input  information  is  available. 
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3.  EMPIRICAL  RELATIONSHIPS 


The  initiation  of  any  hydrometeor  microphysical  modeling  requires  input  informa¬ 
tion  to  quantify  the  total  number  of  particles  and  to  characterize  the  manner  that  the 
particles  are  distributed,  such  as  in  the  form  and  limits  of  the  distribution.  The  3-D 
model  provides  the  liquid/ice  water  "content”  or  density,  M  (Eq.  (8)),  and  ambient 
temperature  (T°C)  of  hydrometeors  contained  within  a  specific  spatial  volume.  M  is 
an  indirect  measure  of  number  density  since  it  is  dependent  upon  both  the  number 
and  sizes  of  the  hydrometeors  being  considered.  Utilization  of  the  hydrometeor  model 
specifies  an  exponential  form  of  particle  distribution;  however,  any  given  M  may  be 
represented  by  a  multitude  of  exponential  distributions  with  various  slopes  and  num¬ 
bers  of  particles. 

With  a  second  measure  of  a  particular  distribution,  such  as  Z,  one  may  gain 
knowledge  of  A  through  the  distributed  values  of  M  and  Z ,  since  M  is  a  function  of  the 
drop  diameters  to  the  third  power  (Eq.  (6))  while  Z  is  a  function  of  the  diameters  to 
the  sixth  power  (Eq.  (11)).  Even  with  this  additional  information,  the  total  number 
of  hydrometeors  may  vary  considerably  depending  upon  the  upper  and  lower  trunca¬ 


tion  limits  (d  and  D  )  of  the  distribution, 
m 


Therefore,  since  the  3-D  model  gives  only  the  M  and  T,  additional  parameter¬ 


ization  (Z,  d  and  Dm)  must  be  specified  or  inferred  from  empirical  relationships  to 


supply  the  necessary  information  to  utilize  the  equation  set  of  Section  2. 

For  any  single  distribution,  M  and  Z  can  be  expressed  as  a  proportionality  fac¬ 
tor  by 


k  =  M  Z 


-0.5 


-3  -1.5 

g  mm  m 


(23) 


This  k  factor  will  vary  depending  upon  the  number  and  size  distribution  of  the 


drops  contained  within  the  diameter  limits  of  d  to  D  .  Although  k  is  determined  from 


m 


calculations  of  M  and  Z  that  are  functions  of  drop  sizes,  the  concept  may  be  utilized 


with  snow /ice  hydrometeors  since  k  is  relatively  insensitive  as  to  the  method  used  to 
deduce  equivalent  melted  diameters  from  the  frozen  particles.10 


♦In  many  cloud  models  that  predict  mixing  ratios  (r),  M  is  often  calculated  as  the 


product  rp  ,  where 


pQ  is  the  dry  air  density. 


10. 


Plank,  V.G.,  Berthel,  R.O.,  and  Barnes,  A. A.,  Jr.  (  1980)  An  improved  method 
for  obtaining  water  content  values  of  ice  hydrometers  from  aircraft  and  radar 


data,  J.  Appl.  Meteor.  19:1293-1299,  AFGL-TR  -  81-0011 ,  AD  A094328. 
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vv*- 


This  relationship  is  valid  for  any  single  hydrometeor  distribution  no  matter  what 
shape  or  form  since  k  can  also  be  written  as 


k  =  C  F  N, 


0.5 


g  mm~3  m 


(24) 


where 


C=f  x  10-® 


(25) 


and  F  is  a  "form  factor"  that  describes  drop  diameter  apportionment  relative  to  each 


other  and  to  the  total  number  concentration  of  the  complete  distribution.  These  ex¬ 
pressions  relating  k  and  N_  have  been  described  and  discussed  elsewhere. 

Id  ^ 
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Past  investigations  have  shown  that  k  varies  with  T  as 


k  =  .0122  e 


-.0474  T 


-3  -1.5 
g  mm  m 


(26) 


for  data  taken  over  Wallops  Island.  VA  and 


k  =  .0143  e 


-.  0291  T 


-3  -1.5 
g  mm  m 


(27) 


for  data  acquired  over  the  Kwajalein  Atolls  in  the  Pacific  for  T  between  -70  and  20°C . 

Equations  (26)  and  (27)  are  derived  from  all  data  available  in  the  data  sets  and 
considerable  scatter  is  evident  since  they  include  individual  aircraft  samples  of  vary¬ 
ing  time  durations  in  different  weather  situations.  These  variabilities  affect  the  k  fac¬ 


tors  as  ( 1)  the  shape  of  the  hydrometeor  distribution  becomes  more  exponential  with 
13 


longer  averaging  time,  and  (2)  storms  experiencing  atmospheric  turbulence  produce 
more  mixing  of  hydrometeor  types  at  comparable  temperature  levels.  No  compensation 


11. 


Plank,  V.G.  (1977)  Hydrometeor  data  and  analytical-theoretical  investigations 
pertaining  to  the  SAMS  Missile  Flights  of  the  1972-73  season  at  Wallops  Island, 
Virginia,  Environmental  Research  Papers,  No.  603,  AFGL/SAMS  Report  No.  5, 
AFGL-TR- 77-0149,  AD  A05U92,  239  pp. 


12. 


Dyer,  R.M.,  Berthel,  R.O.,  and  Izumi,  Y.  (1981)  Techniques  for  measuring 
liquid  water  content  along  a  trajectory.  Environmental  Research  Papers, 
No.  733,  AFGL-TR -81- 0082,  AD  A102922,  52  pp. 


13. 


Plank,  V.G.,  Berthel,  R.O.,  and  L.V.  Delgado  (1980)  The  shape  of  raindrop 
spectra  for  different  situations  and  averaging  periods,  J.  Rech.  Atmos.  14:301- 
309,  AFGL-TR -8 1-0008,  AD  A094877.  ^ 
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can  b«  made  for  the  varying  time  durations  of  the  original  data,  but  the  selection  of 
samples  in  accordance  with  storm  type  serves  to  mitigate  the  effects  of  different  weather 
situations.  Since  the  present  usage  of  the  3-D  model  is  for  convective  storms,  the  k 
data  were  reanalysed  to  form  an  equation  tailored  to  this  type. 

The  Wallops  data,  consisting  mainly  of  stratified  storms  and  cloud  situations,  were 
disregarded .  The  Kwajalein  data,  on  the  other  hand,  are  all  from  convective  cases. 
However,  these  data  include  one  particularly  intense  storm  that  produced  a  k-T  rela¬ 
tionship 


k  =  .0174  e  *042  T  gnm^m  1,5 


that  is  quite  different  from  the  more  normal  situations  (Figure  3).  The  remaining 
storms  gave 


k  =  . 00963  e 


-.031  T 


g  mm  ^  ra 


thus  Eq.  (29)  was  chosen  for  use  in  the  3-D  model  analysis. 

When  M  and  T  are  known,  Eq.  (29)  may  be  substituted  into  Eq.  (23)  and  solved 
for  Z  to  determine  the  equivalent -radar  reflectivity  for  a  particular  hydrometeor  dis¬ 
tribution  as 


.00963  e 


M _ 

-.031  T 


6  -3 

mm  m 


Temperature  may  also  be  associated  with  the  upper  truncation  limit  as  it  can  be 


indirectly  linked  to  a  combination  of  Dm  and  A  ,  as  used  in  Eqs.  (10)  and  (15),  when 


considering  an  exponential  form  for  number  density. 


The  combined  parameter,  D  A  ,  has  previously  been  related  to  hydrometeor 


11  14  15 

type  with  a  value  of  8.2  determined  for  widespread  rain  ,  9  for  large  snow  and 


the  approximate  values  of  10.5  and  12  assigned  to  small  snow  and  ice  crystals  (bullet 


14.  Plank,  V.G.,  and  Berthel,  R.O.  (1982)  A  descriptive  double-truncated  exponential 
model  for  hydrometeors  of  precipitable  size,  Conf.  on  Cloud  Physics,  Chicago. 
Nov  15-18:190-194,  AFGL-TR -82-0347,  AD  A122036. 


15.  Berthel,  R.O.  (1980)  A  method  to  predict  the  parameters  of  a  full  spectral  dis¬ 
tribution  from  instrumentally  truncated  data.  Environmental  Research  Papers. 
No.  689.  AFGL-TR -80- 0001,  AD  A085950,  48  pp. 
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Figure  3.  Aircraft  Acquired  k  Factors  and  Temperatures  from  the  Kwajalein 
Data 


rosettes).  When  the  averaged  T's  from  the  aircraft  measurements  of  the  typical 
Kwajalein  storms  are  associated  with  the  DmA 's  of  the  various  hydrometeor  regions, 
they  generally  conform  to  a  relationship  described  by 


DA  =  8.2  e 

ID 


-.00846  T 


<0°C 


(31) 


with 


DA  =8.2 
m 


>0°C. 


(32) 
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It  is  assumed  that  Eqa.  (31  and  32)  are  valid  up  to  the  raindrop  break-up  sise  of 
~6  mm.  At  the  point  where  D  reaches  6  mm,  A  equals  1.367  mm”*  and  corresponds  to 

03 


a  particular  mass.  Mg .  Since  k  continues  to  decrease  with  increasing  T,  most  likely 
reflecting  evaporation  and  coalescence  processes ,  MQ  is  dependent  upon  T  and  can  be 
calculated  using  Eq.'s  (8),  (10)  and  (29)  (the  value  used  for  d  is  relatively  unimpor¬ 


tant  in  these  calculations  of  large  M  and  Z's).  The  M  at  the  point  where  Dm  reaches 


break-up  sise  forms  a  direct  relationship  with  T  as 


Mb  =  6.144  e 


-.062  T 


As  M  increases  in  excess  of  M_  with  D_  held  constant  at  6  mm,  the  definition  of 

O  09 


D  A  in  Eq.  (31)  and  (32)  is  no  longer  valid  and  is  replaced  by 
m 


D  A  =  6 
m 


in  the  calculations  of  rM  (Eq.  (10))  and  rz  (Eq.  (15)). 


The  Z  and  DfflA  values  determined  from  the  above  relationships  may  now  be  com¬ 


bined  with  a  specified  lower  truncation  limit,  d,  for  application  to  the  equation  set  of 
Section  2.  Calculations  of  A ,  NQ,  NT  and  Dffl  are  shown  plotted  versus  M  in  Figures 
4  through  7  for  various  values  of  T  and  d . 


4.  DERIVED  RELATIONSHIPS 


The  number  and  sizes  of  hydrometeors  at  some  particular  time  and  place  in  a 
storm's  life  cycle  are  important  requisites  for  the  3-D  model.  Using  this  information, 
the  model  may  then  apply  various  physical  processes  to  extrapolste  future  conditions. 

Assuming  that  the  exponential  distribution  equations  may  provide  reasonably 
accurate  determinations,  the  distributed  number  may  be  deduced  from  Eq.  (1)  and 
the  total  number  from  Eq.  (3).  However,  both  equations  require  prior  knowledge  of 
other  distribution  parameters. 

The  sequence  that  would  normally  be  used  to  parameterize  an  exponential  distri¬ 
bution  from  prior  knowledge  of  M  and  T  is  the  determination  of  Z  (Eq.  (30)),  the 
iterative  calculation  of  A  (Eq.  (20)),  the  evaluation  of  Nq  (Eq.  (21  or  22))  and  the 
application  of  these  parameters  to  the  distributed  N  (Eq.  (1))  or  total  N  (Eq.  (3)) 
equations.  This  procedure  produces  very  accurate  values  for  A,  but  at  a  high  com¬ 
putational  expense.  Although  this  circuitous  procedure  is  highly  complicated. 
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especially  the  iterative  process,  it  may  be  easily  solved  by  computer  operation.  How¬ 
ever,  when  a  multitude  of  individual  distributions  are  required,  as  in  the  everchang- 
ing  water  budget  in  a  3-D  model,  much  valuable  computer  time  would  be  consumed. 
Thus,  we  have  investigated  various  means  of  estimating  A  to  eliminate  the  iteration. 

The  most  straight-forward  method  considered  was  the  determination  of  A  by 
assuming  d  =  0  as  a  first-order  estimate  in  the  calculations  of  rM  and  rz  in  Eq.  (20) 
as 

rM  =  1  -  g^--A  [(Da)3  +  3(D  A)2  +  6D  A  +  6  "I  (35) 

M,o  6  L  n>  ra  m  J 

and 

'2,o  *  1  [<D„A>6  *  «V'>5  *  30'I>mA'4  *  “'V1’ 

+  360(D  A)2  +  720(D  A)  +  72ol  .  (36) 

m  m  j 

We  can  then  rewrite  (20)  and  apply  (23)  to  obtain: 

(37) 

Since  the  product  DmA  is  known  froir  Eq.’s  (31),  (32) ,  or  (34)  and  k  is  known  from 
Eq.  (29),  both  as  functions  of  temperature,  evaluation  of  A^  using  this  expression  is 
very  simple. 

We  then  evaluated  the  assumption  d  =  0  by  comparing  values  of  A's  obtained  by 
the  simple  solution  with  those  obtained  by  the  more  accurate  iterative  method. 

When  the  A's  estimated  from  the  simple  method  (called  A£)  are  compared  with 
those  calculated  by  iteration  using  a  specified  nonzero  lower  truncation  limit  (referred 

_3 

to  as  A  ) ,  they  agree  favorably  at  higher  values  of  M  (  >0.1  gm  )  and  also  at  the  high¬ 
er  temperatures  as  shown  in  the  plots  for  d  =  .05  and  .1  mm  in  Figure  8.  The  agree¬ 
ment  becomes  noticeably  better  as  the  specified  value  for  d  becomes  smaller,  and  when 
d  =  0  the  ratio  of  A  g/  A  becomes  unity.  The  discrepancies  between  A  £  and  A  for  d  >  0 
are  greater  at  larger  d  values,  smaller  M's  and  lower  temperatures. 
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The  estimates  provided  by  assuming  d  =  0  are,  in  most  cases,  good  for  situations 
that  may  be  classified  as  being  normal.  For  example,  if  one  were  to  consider  a  situ- 
at  ion  where  M  =  .01  g  m  at  — 40°C,  A  ^  would  equal  28.8  mm-  .  If  under  these  condi¬ 
tions  the  lower  truncation  limit  were  designated  as  being  .  1  mm ,  A  would  be  a  factor 

-l  11 

of  .835  of  the  A  calculated  by  iteration  (34.5  mm  ).  However,  at  -40°C  one  would  ex¬ 
pect  particle  distributions  to  have  lower  values  of  d ,  such  as  .  04  mm ,  where  A  now 
becomes  a  factor  of  .99  of  the  calculated  A  of  29. 1  mm'1.  Thus,  this  method  is  adequate 
for  the  estimation  of  A  when  considering  a  normal  situation.  Non-normal  cases  (i.e. 
cases  with  large  d,  small  M,  and  low  T)  will  either  be  required  to  revert  to  the  trial- 
and-error  solution  or  to  use  to  some  other  correction  procedure. 

Figure  9  shows  the  M  and  T  values  that  define  a  factor  (f)  of  .99  from  Ag/A  for 
various  values  of  d .  Any  combination  of  M  and  T  thaf.  fall  above  the  line  of  designated 
d  will  give  a  A  within  1%  of  the  calculated  value.  Similar  relationships  as  shown  in 
Figure  9  can  be  derived  for  other  factors  and  may  be  expressed  in  terms  of  M ,  T  and  d 
as 


f  = 


(d3  e^066  T) 


-.94 


-  8.26 


-8.22 


(38) 


When  f  ^  .5,  Eq.  (38)  may  thus  be  used  to  adjust  aA£  that  was  determined  by  assum¬ 
ing  d  =  0  by 


A  =  AE/f 


(39) 


to  produce  a  A  that  will  agree  within  1%  to  one  calculated  from  the  equation  set  of  Sec¬ 
tion  2.  Note  that  the  use  of  Eq.  (38)  allows  d  to  be  determined  independently.  For 
example,  smaller  values  of  d  may  be  used  at  altitudes  where  clouds  are  present,  with 
larger  values  being  employed  below  cloud  base  to  represent  the  effects  of  evaporation. 

The  values  of  Nq  may  then  be  found  by  applying  the  corrected  Ag  to  Eqs.  (10) 
and  (21)  and  to  Eqs.  (3)  and  (5)  for  the  determination  of  N,g. 

The  d  =  0  method  for  estimating  A  is  only  applicable  to  distributions  where  D  is 
smaller  than  the  break-up  size  of  rain  drops  (Dm  <  6  mm)  since  Dm  is  a  constant  for 

M  quantities  above  M„  (Eq.  (33)).  However,  estimations  of  A,  when  D  =6  mm,  may 
t>  m 

be  found  by  the  equations 


=  2.293  e 


-.084  M  e 


.062  T 


mm 


-1 


(40) 
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Figure  9.  The  M  and  T  Values  that  Define  a  Factor  (Ag/A)  of  .99  for 
d  =  .02,  .04,  .06,  .08  and  .  1  mm 


for  M  values  less  than 


M  _  . .  -.  062  T  „  —3 
M  =  11  e  g  m 


A  E  =  1.  793  -  .  0802  M  e '  062  T  mm  1 


when  M  is  equal  to  or  larger  than  Eq.  (41).  The  Ag's  determined  from  Eqs.  (40)  and 
(42)  are  within  1%  of  the  calculated  values. 
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5.  REMARKS 


The  concept  of  hydrometeor  modeling  presented  in  this  report  involves  a  blending 
of  theoretical  and  empirical  relationships.  It  should  be  considered  as  a  starting  point 
or  corner  stone  for  the  building  of  realistic  hydrometeor  distribution  predictions  from 
large-scale  meteorological  models  such  as  the  AFGL  3-D  model.  Several  aspects  of  this 
modeling  need  to  be  explored  further. 

Distribution  shape  has  been  investigated  many  times  in  a  variety  of  ways.  The 
large  or  precipitation-sized  drops  generally  show  a  decrease  in  number  with  increasing 
size  and,  in  most  cases,  can  be  reasonably  described  by  an  exponential  curve.  How¬ 
ever,  there  are  distributions  that  do  not  conform  such  as  those  displaying  bi-modal 
tendencies.  Those  non-conformists  can  possibly  be  associated  with  particular  cloud  or 
storm  situations  for  future  modeling  efforts. 

Also,  the  exponential  number  density  is  usually  associated  with  the  precipitation- 

sized  particles.  Since  cloud  hydrometeors  are  heavily  dependent  upon  cloud  type,  the 

shape  of  the  number-density  curve  in  the  small  diameter  range  changes  with  different 
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cloud  situations.  However,  it  has  been  pointed  out  that  a  slight  alteration  by  the 
inclusion  of  Du  and  K  in  the  distribution  function  of  Eq.  (1)  to  form  a  gamma  distribu¬ 
tion  as 


N=KNQDUe'BA  No.  m  3  mm1,  (43) 

where  K  =  1/mm  n  ,  can  conveniently  change  the  curve  shape  especially  in  the  small 
size  range.  Figure  10  shows  the  varying  shapes  that  result  from  different  values  of  u 
in  Eq.  (43)  with  u  =  0  being  an  exponential  distribution.  Again,  association  of  differ¬ 
ent  cloud  and  storm  situations  with  this  type  function  may  result  in  more  realistic  dis¬ 
tribution  predictions. 

The  relationships  of  T  with  DmA  and  the  k  factor  have  been  investigated  for  two 
locations.  Wallops  Island  and  Kwajalein,  and  deserve  further  study.  Associations  of 
hydrometeor  type  with  both  DmA  and  k  have  been  employed  in  the  past  with  varying 
amounts  of  success  but  the  concept  of  relating  these  parameters  with  T  has  only 


16.  Ulbrich,  C.W.  (1984)  The  effects  of  drop  size  distribution  truncation  on  rainfall 
integral  parameters  and  empirical  relations,  J.  Appl.  Meter.  24:580-590. 
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Figure  10.  Changes  in  the  Number  Density  Curve  from  Different 
Values  of  u 
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recently  been  a  serious  consideration.  A  reanalysis  of  past  aircraft  -  acq  uired ,  PMS  1-1 
data17  is  currently  underway  at  AFGL  in  an  attempt  to  verify  or  tighten  these  relation 
ships  and/or  to  derive  equations  for  application  to  different  weather  situations. 


17.  Knollenberg,  R.G.  (1970)  The  optical  array:  an  alternative  to  scattering  or 

extinction  for  airborne  particle  size  determination,  J.  Appl.  Meteor. ,  9(No.  1 
86-103.  ^ 
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